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1. Introduction 

The Dedekind sum can be defined for two relatively prime positive integers 
a, b by 

1 r-i irka irk 
s (a, b) = — > cot — — cot — — . 
v ' 4b ^ b b 
fc=l 

These sums appear in various branches of mathematics: Number Theory, 
Algebraic Geometry, and Topology; they have consequently been studied 
extensively in various contexts. These include the quadratic reciprocity law 
f |13|). random number generators ( |12j). group actions on complex manifolds 
and lattice point problems (|14j. 0). Dedekind was the first to show 
the following reciprocity law (0): 

S (a,6) +S (6,a) = -1 + 1(^ + 1 + ^) (1) 

He was led naturally to this reciprocity law by considering the ^-function 
t/(t) = e^2~IT^ =1 (l — e 27r4m/r ) on the complex upper half plane and trans- 
forming it under the action of the modular group SLi (Z). 

GauB's law of quadratic reciprocity, for example, follows easily from (TjQ) 
(EH; We note that s(a,b) = s(a mod b,b). Combining this with the 

reciprocity law l|T]). one obtains a polynomial-time algorithm for computing 
s(a, b), similar in spirit to the Euclidean algorithm. From this point of 
view, it is not surprising (though not obvious) that s(a, b) can be expressed 
efficiently in terms of the continued fraction expansion of a/6 (jHJ, |19j^. 

Rademacher was one of the pioneers in the use of Dedekind sums (|17|): 
in fact, he found several proofs of (JJ) ( JH])- We present yet another proof, 
which establishes a simple connection with lattice point enumeration in poly- 
topes. The reciprocity law (TjQ) follows readily once the reader is familiar with 
the computation of the coefficients of the Ehrhart polynomial for a lattice 
polytope. 



1 This paper appeared in American Mathematical Monthly 106, no. 5 (1999), 459-462. 



2. Counting Lattice Points 

Let Z n C M n be the n-dimensional integer lattice, and V an n-dimensional 
lattice polytope in W 1 . So V is a compact simplicial complex of pure dimen- 
sion n whose vertices lie on the lattice. For t G N, denote by L(V, t) the num- 
ber of lattice points in the closure of the dilated polytope tV := {tx : x G V}. 
Ehrhart ([5]) proved that L(V, t) is a polynomial in t of degree n. Moreover, 



Here, Vol(dV) denotes the surface area of V normalized with respect to the 
sublattice on each face of V, and xCP) ^ s the Euler characteristic of V. We 
note that, for convex polytopes V, x("P) = 1 

In this paper, we focus on the case M 2 , where Ehrhart's result is known as 
Pick's Theorem [I]): For a convex lattice polytope V G M 2 , 



where A is the area and B the number of boundary lattice points of V . 

In the general case, the other coefficients of L(V, t) are not as easily acces- 
sible. In fact, until quite recently a method of computing these coefficients 
was unknown. There has been recent progress in this direction (pQ, 0, 
|10| , ) ; Diaz and Robins ( [H] ) found a way of proving a cotangent repre- 
sentation for the generating function YltLo ^ / (^ > ' t)e~ 27TSt , thereby deriving a 
formula for the Ehrhart coefficients of L(V, t). For our purposes, the follow- 
ing result (a straightforward consequence of Corollary 1 in 5 ) is sufficient: 

Theorem. Let V denote the simplex in W 1 with the vertices 
(0, . . . ,0), (ai,0, . . . 0) , (0, a 2 , 0, . . . ,0) , . . . , (0, . . . ,0, a n ), where oi, . . . , a n G 
N are pairwise coprime. Denote the corresponding Ehrhart polynomial by 
L{V,t) = E]= Cjt j . Then c m is the coefficient of s ( m + 1 ) in the Laurent 
expansion at s = of 






where p = a\ ■ . . . ■ a. 
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Note that the appearance of cotangent products in this result leads us to 
expect Dedekind sums in some form within the coefficients of the Ehrhart 
polynomial, thus also within the formulas for the number of lattice points 
in simplices. In fact, the nontrivial cases of dimension three (|15|) and four 
([15]) involve classical Dedekind sums. Both formulas can easily be obtained 
through the above Theorem. 

We will use this result in an indirect way. Precisely, we will compute Co 
according to the above theorem, and make use of the fact that co = xi'P) = 
1. Dedekind's reciprocity law Q will follow from this idea if we consider 
the case of dimension n = 2. 



3. Proof of the Reciprocity Law 

According to the theorem, we have to find the coefficient of s" 1 of the 
Laurent series at s = of (note that (a, b) = 1) 



ab 

7T 



— - ( 1 + coth — (s + irj] (1 + coth — (s + irj] ( 1 + coth — (s + ir 
Aab^K a J \ b K J \ ab x 

r=l 

(2) 

The Laurent expansion of each factor depends on r: 

i, j.l ^ / , ■ \ / S c :=^s- 1 + l + §- c s + 0(s s ) iic\r 

1 + COtll — ( S + IT ) = < „ T i iriv s~\ / \ -r 1/ 

c y ' \ R c ■= 1 + coth ?f + 0(s) if c /r 

To keep track of the various cases, we introduce the notation 

J 1 if c|r 

Xc= { if c|r, 

so that we can write 1 + coth ^(s + ir) = S c x c + Rc (1 — Xc), and © becomes 

ab 

^2 ( S aX°- + R aO-- Xa)) (S b Xb + -R& (1 - Xb)) (S a bXab + Rab (1 ~ Xab)) ■ 
r=l 



Now, expand this into all 8 terms, and consider each summand according to 
the number of S c factors: 
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1. Terms with one S c factor are 

SaXaRb (1 - Xb) Rab (1 ~ Xab) = S a R b R ab Xa (1 ~ Xb ~ Xab + Xab) 

= S a R b R a b {Xa ~ Xab) (3) 

and, similarly, 

Ra (1 - Xa) SbXbRab (1 ~ Xab) = R a SbR a b (Xb ~ Xab) ■ (4) 

The summand with S ab is zero (note that XaXab = XbXab = Xab-, an d 
XaXb = Xab)- To compute the contribution of (J3J), note that the sup- 
port of Xa — Xab in {1, ... , ah] is {ka : 1 < k < b — 1}; thus its contri- 
bution to © is 



a sr^ ( , 7rika\ ( , nika\ 

- > 1 + coth — — 1 + coth — — 

vr ^ V & / V ab J 

1 



7rA;a\ / irk 
i cot — — 1 — i cot — - 



7T 

Aal> 



46 ^ I 6 7 V 6 

fe=l v 7 v 

1 ^-i 7rfca irk 

= — > 1 - cot — — cot — + % . . . 
46 ^ 6 6 

fc=i 

= 4-46- 5(a ' 6) - 

The imaginary part in the preceding sum has to be zero, because the 
original generating function is real. Similarly, Q gives a contribution 
of I - £ ~s(b,a). 

2. There are no terms with two S c factors, because 

SaXaSbXbRab (1 ~ Xab) = S a S b R ab Xab (1 ~ Xab) = 

and 

SaXaRb (1 - Xb) SabXab = S a R b S a bXab (1 ~ Xb) = 0. 

3. Finally, the term S a x a SbXbS a bXab = S a S b S a bXab has support {ab}, and 
gives a contribution of 

7T / a b 7T 
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Adding all contributions, we arrive at 

l = co = i + -^- + - + -j- S (a,6)- 5 (6,a), 
the desired reciprocity law 

The same method applied to dimension n = 3 does not give any further 
results. However, for n = 4, higher dimensional Dedekind sums ([201) appear 
within the computations, so that this case is likely to provide new results. I 
will discuss this matter in a future paper. 
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